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ABSTRACT 

A single phenomenological constitutive equation is derived theoretically from 
first principles and applied to aluminum, tin and lead. The theory is based on 
deformation kinetics of steady creep in which the fundamental mechanism is atomic 
transport over potential barriers whose conformation is distorted by the applica- 
tion of a stress field. 

The form of the functional dependence of barrier distortion and stress over 
the entire temperature range is found to be a sigmoidal curve which tends to 
straight lines of a unit slope in the small and high stress regions. With this 
form of barrier distortion, the constitutive equation can predict very well the 
steady creep behavior of aluminum, tin and lead over a wide range of temperature 
and stress. 

1. Introduction 


Experimental results on high temperature creep of pure metals and solid solution 

alloys during past decades, [1,2,3] fall into two main categories: those establishing 

a relationship between steady-state creep rate and stress under constant temperature 

on one hand and a relationship between constant creep rate and temperature under 

constant stress, on the other. The functional dependence of the constant creep rate 

e on the stress a under constant absolute temperature T may be divided basically into 
s 

three regions whose boundaries depend on the material itself. In the low stress 
region, e s is almost linearly proportional to a. Hence it is called a newtonian- 
like viscous flow region. In the intermediate stress region, e appears to be 
proportional to J 1 , where n is a temperature dependent material parameter. The value 
of n lies predominatly between 4 and 7 for pure metals ^nd between 3 and 5 for solid 
solution alloys. In the high stress region, e s is proportional to the exponential 
function of 8a, here 8 is a temperature dependent material parameter. The functional 
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dependence of e on T under constant o is assumed to be governed by a type of 
Arrhenius relation. As a result the activation energy of creep can be found by a 
temperature cycling technique or the slope of the line in the Arrhenius plot of 
log^ vs T \ The values of activation energy thus found are very close to those 
of self diffusion in pure metals or the diffusion of one of the predominant elements 
in solid solution alloys. However, the slope of the Arrhenius plot is, in general, a 
function of stress and temperature. In addition the activation energy calculated 
from temperature cycling technique is, in general, a function of stress and strain. 

Various theories of creep have been proposed in. recent times. These fall 
basically into two broad categories: phenomenological and micromechanical , the final 

aim being, of course, a macroscopic constitutive equation. The first category 
includes theories that are strictly empirical [2-4], others that are mathematical [5], 
and others still which are "quasi" physical such as the internal state variable 
theories, a typical example being reference 6. In the second category fall theories 
in which the underlying micromechanisms are vacancy diffusion, dislocation climb and 
microcreep [7]. In the latter category belong also the absolute reaction rate theory 
by Eyring [8] and the very recent deformation kinetics theory of creep by the 
authors [9]. 

Micromechanical theories, where vacancies or dislocation are the building blocks, 
need more than one mechanism to describe the experimental phenomena over a wide range 
of temperature and stress. While for practical purposes this is not a disadvantage, 
one wonders if a single appropriate atomic mechanism cannot be found which describes 
steady creep phenomena over the entire range of stress and temperature. 

In reference 9, we found reason to believe that this might be possible. One 
single constitutive equation was shown to predict very well the steady creep behavior 
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of AISI 316 stainless steel, pure polycrystalline aluminum and copper over a wide 
range of temperature and stress, above about lOO^cm^. The vehicle for this specific 
constitutive equation is deformation kinetics. The fundamental mechanism is the 
transport (diffusion) of atoms over potential barriers whose conformation is distorted 
by the application of a stress field. 

In the case of one-dimensional flow, of interest here, the central element of the 
theory is the relation between the barrier distortion m and the free energy gradient 


_ _ where q (an internal variable) is the statistical average of the displacement of 
9q 

atoms in motion facing a specific barrier . In the application of the theory [9] to 

uniaxial stress fields where the stress was above circa lOO^^m^ a linear relation 

between cj and - — sufficed but proved inappropriate for lower stress levels, 
dq 

Evidently the task at hand is to find an appropriate relation that applies to all 
stress levels but the form of the relation does not negate the fact that we are 
dealing with a single mechanism of atomic diffusion over energy barriers. This is 
done in Section 3. 

2. Brief Review of the Theory 
Particle Equations 

Let N be the number of particles whose motion is impeded by a barrier of height 
£ q and to the distortion of the barrier due to the application of the stress field. 

See reference 9 for details. Then the number of particles N ; partaking in the net 
motion is given by equation (2.1) 

N' = 2Nexp(-e o /kT)sinhWkT) (2.1) 

where k is the Boltzmann constant and T the absolute temperature. Assuming a "square 
sinusoidal" barrier shape the average time x taken by the atoms to climb the barrier 
is given by equation (2.2) 
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( 2 . 2 ) 


H a . -1, , 

t = —p=~ F (2, sm k ) 
Tr/e o 


where a is the barrier width, F is the complete elliptic integral and 


■ 


(2.3) 


Rate (Evolution) Equation for a 

The average velocity q of the atoms crossing the barrier is given by equation 

(2.4) 



(2.4) 


Use of equations (2.1), (2.2) and (2.4) gives the desired relation between the 
average velocity and the barrier distortion: 

Jr. [2e 

q = exp (-£ / kT) sinh (to/kT) (2.5) 

r o 


A convenient representation for F is the following 

F = \ log (16e 0 /u>) (2.6) 

For o<o)/e o <0.35 the maximum error is less than 5% [9]. If there exists n barriers to 
the motion, each of height s o with distortion then equation (2.5) applies to each 
such barrier. 

However in steady creep only the highest barriers come into play, the lowest ones 
having already been climbed by the atoms in the course of the deformation. Thus one 
internal variable suffices (approximately) to represent the effect of these barriers, 
if, indeed more than one is actual active. Otherwise the representation is exact. 
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Free Energy Representation 


In general 

ip = i|i( e .q,T) (2.7) 

where in the present work T is constant. The free energy is the potential energy 
stored by virtue of atoms being displaced within potential wells. The mean 
displacement generated as a result is directly related to the elastic strain. For 
instance in the case of a unidirectional equispaced atoms the elastic strain is 
exactly equal to the atomic displacement divided by the lattice spacing. Assuming 
parabolic wells, the potential energy is proportional to a quadratic function of the 
displacement, leading to the conclusion that the free energy is a quadratic function 
of the elastic strain [9]. 

To relate the above discussion to equation (2.7) we write \p in the quadratic form 

= h A n e 2 + A 12 eq + h A 22 q 2 (2.8) 

and insist that it is a perfect square, so that the squared linear term can then be 

2 

identified with the elastic strain. This is possible if A^ = A ^22' T ^ e P r: *- nc: *-Pl e 

of thermodynamic stability requires that be positive definite. This implies 

2 

An > ^22 >0 ’ ^12 -^ 22 ^ 11 < °* However, the last inequality can be relaxed and set 

into equality for the purpose of steady state creep in which the metal exhibits a 
fluid equilibrium configuration [9]. As a restilt, equation (2.8) becomes 

^ = %A( e - Bq) 2 (2.9) 

where A = A ^ and B = -A^/A. Thus since equation (2.9) is the mathematization of the 
statement at the end of the last paragraph e -Bq must be identified as an elastic 
strain. Note that A and B may be and are, in general, functions of temperature. 

To obtain the desired analytical expression for creep we appeal to a 
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fundamental relation of irreversible thermodynamics according to which the stress is 


the gradient of free energy with respect to the strain i.e., o = dip/ de. Thus 
a = A(e- Bq) . Furthermore, as a result of equation (2.9), 9tp/3q --Bo. Since 
during creep the stress is constant, it follows that e = B q. The strain rate can 
then be obtained from equations (2.5) and (2.6), i.e., 

-£ /kT 


£ = 2/2e irBe 
o 


sinh(m/kT) 


( 2 . 10 ) 


where to is now a function of a and T. 

In the next section we will use equation (2.10) to predict the steady creep of 
aluminum, tin and lead, particularly under very high temperature and low stress. 
However before this can be done the relation between the internal force - 3^/ 3q and 
to must be established. As noted above - d^/dq - Bo. The problem is therefore reduced 
to finding the relation between o and to, in this particular case. 

3. Application of the Theory to Aluminum 

Let the relation cr (to) between a and to or conversely, w(cr) be known. Specifically 

let 


a = 2 (w) ; to = 0 (a) (3.1a,b) 

Substitution of equation (3.1b) in equation (2.10) gives a constitutive relation 

-e /kT 


e = e 


J (a,T) 


(3.2) 


In an inverse fashion, given the experimental relation between a and e at constant T, 
one may then use equation (2.10) to deduce the relation between a and to, i.e., the 
function 2 (to) . 

The function of 2 (to) in the case of aluminum can be found from figures 1 and 2. 

It may be seen that for higher values of stress the relation of a and log n £ is 

XU s 
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linear. Thus, in this range, the linear relation implies that 


= kTK, 


i 1 - 


lcrl 


(3.3) 


where a 1 is a threshold stress below which equation (3.3) is not applicable, and K~ 
o 

is the slope of the straight portion of the curve. Equation (3.3) was the basis of 
the study in reference 9, where it was shown that the linear relation persists over 
a wide range of temperature with the proviso that and aj are temperature 
dependent. In this case^(a,T) has the form 


J*(0, T) = 


K° sinh K„( a- a* ) 
1 2 o 


log 


16 e 

-kx 2 - / K 2 C- °o > 


(3.4) 


where K° = 2^2e ^B. The determination of the constants K° , K 0 and o was discussed 

1 o 1 2 o 

at length in reference 9. Note that equation (3.3) implies that the distortion is 
linearly related to the internal force, i.e., 

to = kTC (Q-Q°) (3.5) 

where Q is defined as - S^/Sq and thus equal to Ba, Q° = BaJ and C is the 

coefficient of proportionality. 

As shown in figures 1 and 2, equations (3.2) and (3.4) predict quite well the 

experimental data in references 10 and 11, except the temperature at 920°k. The 

corresponding values of K°, K^, C and aj are shown in figures 3 and 4; e q = 34 

KCal/mole which is the value of activation energy of self diffusion. It is seen 

that serious deviations begin to arise below a stress level of about a Of 

o 

course this is to be expected in view of equation (3.3). Below , equation (3.3) 
no longer holds. 

To determine £2(a) is this region we recall equation (2.10), which because of 
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the smallness of w/kT we can write in the approximate form 


= K? e 


-e 0 /kT 


ut/kT 

log(16e Q /a)) 


(3.6) 


and note that insofar as this region is concerned the experimental data at 920°k 


indicate a linear relation between log^Q £ s an d log-^g ~ ° Q ) » w ^ en th e net stress 

It *1 

a - o q ) , is very small [11,12]. As indicated in reference 11, o q = 3psi , below 
which creep was not measurable. Since log log (16e o /w) is an insensitive function of 
m in this small net stress region, the above observation suggests the following 
relation between u) and a: 


II 

log 1Q 6m = lo S 10 A 0 + lo S 10 (a " a 0 ^ (3-7) 

If 

where A and a are at most functions of temperature; R = 1/kT. These parameters 
o o 

were determined respectively from the intercept and the constraint of a unit slope 

-5 H 

of the curve. In the present case A is a constant (4.6x10 ) and o is a 

o o 

decreasing function of temperature (see figure 4) . Indications are that as the 
temperature approaches the melting point (T goes to zero at which point the 

metal exhibits a truly newtonian behavior. 

The form of the functional dependence of w on a over the entire temperature 
range is shown in figure 5. The relationship is sigmoidal tending to a linear form 
in the small and high stress regions (see equations (3.7) and (3.3)). The 
theoretical predictions based on figure 5 are shown in figure 2. 

4. Application of the Theory to Tin and Lead 

The procedure of Section 3 is repeated here without change. It has been shown 
by a three-dimensional argument [13] that the constitutive equation (2.10) applies 
to pure shear without change in form. A comparison between theory and experiment 
is shown in figures 6 and 7. A further treatment will be the subject of a more 
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extensive article to appear at a later date. 


5. Conclusions 

In this paper we apply the theory of deformation kinetics to aluminum, tin and 
lead and show that equation (2.10) suffices to predict accurately steady creep 
behavior over a wide range of stress and temperature. We may conclude that one 
micromechanism, that of atomic transport over potential barriers whose conformation 
is distorted by the application of a stress field, is sufficient to account for the 
steady creep process in the entire range of temperature and stress. 

References 


1. Dorn, J.E., "Some Fundamental Experiments on High Temperature Creep," J. Mech. 
Phy. of Solids, 3, 1954, 85-116. 

2. Garofalo, F., "Fundamentals of Creep and Creep-Rupture in Metals," MacMillan Co. 
New York, 1965. 

3. Sherby, O.D. and Burke, P.M., "Mechanical Behavior of Crystalline Solids at 
Elevated Temperature," Progress in Material Science, ed. B. Chalmers and Hume- 
Rothery, Vol. 13, 1968, 325-390. 

4. Nadai, A. and McVetty, P.G., "Hyperbolic Sine Chart for Estimating Working 
Stresses of Alloys at Elevated Temperature," Proc. ASTM, 43, 1943, 735-748. 

5. Findley, W.N., Lai, J.S. and Onaran, K. , "Creep and Relaxation of Nonlinear 
Viscoelastic Materials," North-Holland , New York, 1976. 

6. Hart, E.W., "Constitutive Relations for the Nonelastic Deformation of Metals," 
Trans. ASME J. Engng Mat. Tech., 3, 1976, 193-202. 

7. Ashby, M.F., "A First Report on Deformation-Mechanism Maps," Acta Metallurgica, 
20, 1972, 887-897. 

8. Krausz, A.S. and Eyring, H. , "Deformation Kinetics," Wiley, New York, 1975. 

9. Valanis , K.C. and Lee, C.F., "Deformation Kinetics of Steady-State Creep in 
Metals," Inti. J. Solids and Structures, 17, 1981, 589-604. 

10. Servi, I.S. and Grant, N.J., "Creep and Stress Rupture Behavior of Aluminum as a 
Function of Purity," Trans. AIME. 191, 1951, 909-922. 


147 



11. Harper, J.G., Shepard, L.A. and Dorn, J.E., "Creep of Aluminum under Extremely 
Small Stresses," Acta Met. 6, 1958, 509-518. 

12. Mohamed, F.A. , Murry, K.L. and Morris, J.W. Jr., "Harper-Dorn Creep in Al, Pb, 
and Sn," Met. Trans. A, 1973, 935-9A0. 

13. Valanis, K.C., "Deformation Kinetics in Three Dimensions," submitted to 
the Int. J. of Solid and Structures. 


148 



PURE POLYCRYSTALLINE ALUMINUM 



Fig. 2 
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